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Abstract. We discuss a role of the localized tt orbital, which exists around the defect, on 
the defect-induced Kondo effect in graphene by a numerical renormalization group study. 
We find that the localized tt orbital assists this Kondo effect, and the Kondo temperatnre 
is sensitive to the broadening of the localized tt orbital. Secondly, we focus on the negative 
magnetoresistance of this Kondo effect. In the experimental result, it has been shown that the 
negative magnetoresistance is ten times larger than the nsual Kondo effect. In order to clarify 
the mechanism of the ’’magnetic sensitive” Kondo effect, as a first step, we study an orbital 
magnetic field dependence of the localized tt orbital by a tight-binding model with a Peierls 
phase. We find that as the magnetic field increases, the spectral width of the localized tt orbital 
increases and the local DOS at the Fermi level decreases. Since the Kondo temperature is 
strongly dependent of the broadening of the localized tt orbital, it is expected that this Kondo 
effect is sensitive to the orbital magnetic field as observed in the experiment. 


1. Introduction 

After the fabrication of graphene, the Kondo effect in graphene has attracted much theoretical 
attention, because exotic Kondo effects and associated phenomena due to its Dirac conduction 
electrons are expected, such as multichannel Kondo effects, gate-tunable Kondo effects, and an 
impurity quantum phase transition between an unscreened and a screened localized moment state 
[1-6]. As candidates of the realization of the Kondo effect in graphene, the graphene with magnetic 
impurities such as Fe or Co has been investigated by DFT calculations [2] . The possibility of the 
Kondo effect in these systems has been argued by examining the lattice positions, the spin state of 
the impurity atom, and the strength and symmetry of the hybridization between the impurity 
and the conduction electrons. Experiments on graphene with such magnetic impurities have 
been attempted, and spectral signatures of the Kondo screening have been observed in scanning 
tunneling microscope (STM) experiments [?]• In such systems, however, Kondo effect had not 
been observed in resistivity measurements. 

In 2011, the Kondo effect was reported in the resistivity of graphene with defects [8j. 
The defects were introduced by the irradiation of He^ ions with the energy of 500 eV. The 
temperature dependence of resistivity showed logarithmic increase with decreasing temperature 
and saturating behavior at low temperatures. The temperature dependence can be fitted well 
by the universal temperature dependence of the Kondo effect. In the ion-irradiated graphene, 
the negative magnetoresistance was also observed [8]. Magnetic field is applied perpendicular 
to the graphene plane. The contribution of the weak localization is suppressed by the magnetic 
field i? ~ 1 T, and a magnetoresistance due to another mechanism occurs at higher magnetic 


fields. In the case of the Kondo effect, the Zeeman effect of magnetic field can cause such a 
negative magnetoresistance. Prom these observations, the Kondo effect due to magnetic origin 
was concluded in this system. To clarify the origin of the Kondo effect, many theoretical studies 
have been performed [9-15].On the other hand, it also has been proposed that the electron- 
electron interaction in the presence of the disorder is responsible for the logarithmic increase of 
the resistivity at low temperatures [muz]. Although there is the other possible theory, in this 
paper, we focus on the defect induced Kondo effect. 

Until the discovery of the defect induced Kondo effect, the defects in graphene or graphite 
have been studied as a source of magnetism [TB]. In graphite, it has been observed that magnetic 
moments are induced by irradiation of high-energy (2.25 MeV) protons [TU]. In the graphene 
irradiated with ions, the magnetic moments have also been observed [20]. It has been indicated 
that the defect individually behaves as a localized magnetic moment. Such defects have been 
studied in STM experiments ED [22]. These experiments show that the introduced defects 
are point defects (or point vacancies). The defects in graphene have also been investigated 
with transmission electron microscopes (TEMs) [231 ) s-iid reconstruction and deformation of the 
defect are observed. Such a point defect in graphene has been studied with the DFT calculations. 
First, the lattice structure was investigated [2l|- Two possible structures were considered: one is 
a threefold symmetric structure that preserves the original symmetry, and the other is a twofold 
symmetric structure in which one of the C-C distances is shortened. It was clarified that the 
latter is more stable owning to the Jahn-Teller effect, and it has an out-of-plane displacement, 
while the former has a planar structure. Second, the spin state was investigated by using spin- 
polarized DFT calculations [251126j . It was shown that the defect is magnetic and the dangling 
sp^ bond which appears at the defect is responsible for the magnetic moment. 

In graphene, effects of defect scattering on the vr electrons have been studied theoretically, 
and it has been suggested that the point defect induces a localized state of vr electrons around 
the defect whose energy eigenvalue coincides with the energy of the Dirac points [26-29].Such 
a localized vr orbital has been observed experimentally with the STM. In scanning tunneling 
spectroscopy (STS) experiments, the localized vr orbital exhibits a sharp peak in local density 
of states (LDOS) near the defect. 

Although the origin of the Kondo effect in graphene with defects seems to be due to the 
localized magnetic moments on the sp^ orbital at the defects, the role of localized vr orbital 
has not been understood in detail. Recently, although a model where the localized state of vr 
electrons is rather artificially added at the defect [iniiis] has been proposed, more consideration 
on the derivation of the model and the analysis are necessary. In addition, in this defect induced 
Kondo effect, the anomalous negative magnetoresistance has been observed: In the usual Kondo 
effect, the energy scale of a characteristic magnetic field of the negative magnetoresistance is 
the same order as that of the Kondo temperature. On the other hand, experimentally, the 
characteristic magnetic field seems to be one tenth smaller than the value expected from the 
Kondo temperature, which means that this Kondo effect is sensitive to the magnetic field. 
However, this mechanism has not been understood yet. 

Thus, in this paper, we discuss the role of localized vr orbital on the defect-induced Kondo 
effect in graphene on the basis of a numerical renormalization group (NRG) study ED. As a 
result, we find that the localized vr orbital assists the Kondo effect, and the Kondo temperature is 
sensitive to the broadening of the localized vr orbital. Next, in order to clarify the mechanism of 
the magnetic sensitive Kondo effect, as a first step, we study an orbital magnetic field dependence 
of the localized vr orbital by a tight-binding model with a Peierls phase. We show that the LDOS 
on the localized vr orbital decreases as the magnetic field increases. 

In §.2, we develop an effective model with the localized vr orbital and transforms the effective 
model into a one-dimensional representation of the conduction electron states to apply NRG 
method. In §.3, we discuss the effect of the localized vr orbital by NRG study. Finally, in §.4, 


we discuss an orbital magnetic field dependence of the localized state of vr electrons by the tight 
binding model with the Peierls phase. 


2. Localized Moment and its Hybridization with Conduction Electrons 
2.1. Cluster Model for the Defect 

Here, considering a cluster for the defect, we introduce a model for the system. Around the 
point defect, there are three sp^ orbitals (pi which do not form covalent bond, where i indicates 
one of the three sites (i = 1, 2, or 3) around the defect. There are also tt orbitals on carbon 
atoms. 

First, we consider an isolated cluster that consists of three sp'^ orbitals, as a model for the 
localized state of sp^ orbitals around the defect. A tight-binding approximation is applied to 
this cluster. Two of the three transfer integrals between the sp'^ orbitals are —h, and the other 
is —h' {h,h' > 0). The threefold symmetric structure corresponds to the case with h = h', and 
the twofold symmetric one to h' > h, where it is assumed that the distance between sites 2 and 
3 is shorter than the others. 

The energy eigenvalues and the corresponding wave functions of this cluster model are 
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which are depicted in Fig. [TJ There are three electrons on these sp^ orbitals. Since the number 
of electrons is odd, there is an unfilled level and nonzero spin appears. Interpreting the Kohn- 
Sham eigenvalues of the DFT calculations [21] in terms of these tight-binding energy levels, we 
semiquantitatively estimate h'/h to be h'/h ~ 5. When h'/h ~ 5, C± becomes (7+ ~ 1, C_ <C 1, 
and the wave function of the lowest level ifi becomes a bonding state ifi = {(p 2 + while 

the second state becomes 'if 2 = Pi- Thus, pi becomes the highest occupied molecular orbital of 
this cluster, and causes a localized magnetic moment. It should be noted that in a divacancy 
case, four electrons are present in the sp'^ orbitals and the localized spin does not appear. 

Second, we consider the hybridization between this localized magnetic moment on the sp'^ 
orbital and the conduction electrons which consist of tt orbitals. It depends on the lattice 
structure of the defect. When the defect has the planar structure, there is no hybridization 
between the sp^ and vr orbitals, which is schematically shown in Fig. [2]^a). However, with an 
out-of-plane displacement, finite hybridization between the sp'^ orbital pi and the two vr orbitals 
at the other sites j i appears as in Fig. [2]^b). As mentioned above, from the DFT calculations, 
the twofold symmetric structure with an out-of-plane displacement is suggested to be more 
stable owing to the Jahn-Teller effect |24| . 

In the following we consider the structure shown in Fig. [3l where only one sp^ orbital (v = 1) 
is active and it hybridizes with the vr orbitals at sites 2 and 3. (The bonding state of sp^ orbitals 
at sites 2 and 3 is neglected.) We assume that the hybridization term is given by 
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Figure 1. Energy levels and the corresponding wave functions of a cluster model for localized 
states of sp^ orbitals around a point defect in graphene. 



Figure 2. Schematic picture of sp^ (blue) and tt (red) orbitals on the different carbon atoms at 
ri and r 2 . (a) Planar structure, (b) In the presence of out-of-plane displacement of the carbon 
atom. 


where Ojo- is an annihilation operator of a vr electron at site i with spin a =t)i and is an 
annihilation operator of an electron on the sp^ orbital at site i = 1. F is the amplitude of 
hybridization. Taking account of the localized nature of the active sp^ orbital in the defect, we 
assume that the effective Hamiltonian for this orbital becomes 

^def = - ^J^)dlda + Ud\d^d^^di, (4) 

(7 

where e^p 2 , p, and U are the energy level, chemical potential, and Coulomb interaction on the 
sp^ orbital, respectively. The Coulomb interaction term is essential for describing the Kondo 
effect. 

The conduction electron states are formed by vr orbitals in graphene, which are described by 
the following tight-binding model with a nearest-neighbor transfer integral t{> 0), 

^^gra = + h.c.) - pN, (5) 

(bV 

where bj^ is an annihilation operator of a vr electron at site j of the sublattice B, and N is 
the total number of vr electrons. In eq. ([3]), only the vr orbitals of the sublattice A, ai^, appear 
since we have assumed the missing carbon atom to be on the sublattice B. The prime in eq. ([5]) 



















Figure 3. Schematic picture of the electronic state around the defect. Two of the three sp^ 
orbitals {i = 2, 3) form a covalent bond. The circle at the site 1 shows the active sp^ orbital 
and those at the sites 2 and 3 the vr orbitals. V is the amplitude of the hybridization between 
the active sp^ orbital and the two vr orbitals. 


denotes the presence of the point defect. The point defect in tt electrons is represented by 
the zero transfer integrals between the missing site and the neighboring site. Here, we ignore 
the Coulomb interaction in the vr orbitals. Although band structure may be modified by the 
Coulomb interaction in the bulk, the Dirac dispersion is preserved and we assume that eq. (© 
represents such renormalized Dirac fermions. Finally, the effective Hamiltonian of graphene with 
the point defect becomes H = ffgra + -f^def + ^^hyb- 

2.2. One-Dimensional Representation of the Conduction Electron States 

To clarify the conduction electron states which are coupled with the sp'^ orbital at the defect, 
and to apply the NRG method, we transform this effective Hamiltonian Q, ([H, and m into 
a one-dimensional representation. In the usual Kondo problems, a magnetic impurity interacts 
with conduction electrons at the same site. In this case, the conduction electron states are 
expanded by the spherical waves, and states with particular symmetries such as s, p, or d wave 
symmetries are coupled with the impurity. For example, when the interaction between the 
impurity and the conduction electrons is spherical, only the states with s wave symmetry is 
coupled with the impurity. Such states are labeled by the wave vectors in the radial direction. 
As a result, a one-dimensional representation is obtained. However, in the present model, there 
are two centers: the orbital at the defect interacts with the vr orbitals of sites 2 and 3. Such 
a geometry has been treated in the two-impurity Kondo problem [32]. We follow the method 
developed for it [32l |33l [M] , generalizing the method to include the effects of defect scattering 
on the TT conduction electrons. 

To express the hybridization term ([3|) in a one-dimensional representation, we introduce the 
following operator for vr electrons [SS] (Mj ; 

4ia = (6) 

k 

where Sk and \k) are the energy eigenvalues and the wave functions of the Hamiltonian m, 
and is the corresponding creation operator with spin a. \Ai) represents the vr electron state 
localized at the site i{= 2 or 3) of sublattice A. 5(e) is the Dirac delta function. These operators 
are related to the original one Ojo- by 

dea^iQ- — aifj. (7) 

The anti-commutation relations of these operators are calculated as j Oeio-, Og/jg. r = 
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0 , where gij^e) is introduced as 
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The diagonal element guie) corresponds to the LDOS at site i of sublattice A. Owning to 
the twofold symmetry in the present system, the LDOS at the site i = 2 and 3 are identical: 
522 (e) = 533 (e). The off-diagonal elements satisfy 523 (e) = 532 (e), since the wave functions can 
be chosen to be real. Then, from as 2 a and a^sa, the following operators which satisfy the usual 
anti-commutation relations are obtained; 


Cetre — i /o (®£2 (t T deSa) i 

{£) 

C-sao — TT^ (®£2(t O-eScr) ; 

(e) 

where e and o denote even and odd, respectively, and 

A^e (e) = 522 (e) 523 (e) + 532 (e) -b 533 (e), 
A^o (e) = 522 (e) - 523 (e) - 532 (e) -b 533 (e). 

With these operators, the hybridization term ([3]) is expressed as 

^^hyb = J d£N^/‘^{£)(Vcl^^d„ + h.c}j . 
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Since Nf,{e) is essentially the LDOS of vr electrons at the defect, this expression shows the 
importance of the LDOS at the defect, not the total DOS. Note that only one channel of 
the conduction electron Ce^e hybridizes with the sp'^ orbital. This is because the sp'^ orbital 
hybridizes equally with the vr orbitals at sites 2 and 3 as in Fig. [3j When the hybridization is 
not equal, both even and odd states are coupled to the sp'^ orbital. 

Using eqs. (0), dSD, and Q, it can be shown that Ce^e satisfies the commutation relation 
= (e — p,)clae- Thus, neglecting the other channels of the conduction electrons, we 


H, 
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assume that the conduction electron term ([5)) can be written as [31-33] 
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Since decoupled channels of conduction electrons are neglected, this model can be used only 
when the defect contribution is concerned. The transformation above has also been justified in 
ref. [ 2 S] by integrating out the conduction electron degrees of freedom in the functional integral 
representation. 


3. Role of the Localized State of vr Electrons at the Defect: Numerical 
Renormalization Group Study 

In this section, we discuss the effect of the localized state of vr electrons on the Kondo effect. 
The localized state of vr electrons causes a peak in the LDOS on the defect as a function of 
energy at the charge neutrality point. To take into account this peak of the LDOS, we assume 
a model shown in Fig. 01 
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Figure 4. Local density of states of conduction electrons in the model in ^ as a function of 
energy e, where D is a cut-off energy. A peak due to a localized state of vr electrons is included, 
where rc is a width of the peak. 


(a) 
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Figure 5. (a) Temperature dependence of entropy of the electron on the defect sp^ orbital, 
S'def/log2, in the presence of a peak in the local density of states of tt electrons with the width 
w/D = 0.01 at the chemical potential p, = 0, for the hybridization V/D = 0.01, 0.05, and 
0.1. (b) Kondo temperature Tk/D as a function of hybridization V/D for several widths w/D. 
[Reprinted from Taro Kanao, Doctor Thesis, Department of physics. Graduate school of science, 
The university of Tokyo, Tokyo (2014). ] 


where A and w are the intensity and the width of the peak, respectively. The second term 
becomes a Dirac delta function in the limit re —?■ 0. In the case of low density of defects, the 
LDOS at the defect is a Dirac delta function with zero width. On the other hand, in the case 
of finite concentration of defects, the LDOS can have a finite width owning to overlaps of the 
wave functions of the localized states on different defects [381 [39] . Assuming finite widths of the 
LDOS, we analyze this model by the NRG method. The p = 0 case is studied, where the effect 
of the peak of the LDOS is most significant. 

Figure EKa) shows a temperature dependence of entropy of the electron on the defect sp'^ 
orbital, 5def/fe_B log 2, for V/D = 0.01, 0.05, and 0.1 in the presence of finite peak width 
w/D = 0.01. In this NRG calculation, the parameter of the logarithmic intervals. A, is set 
to be A = 2, and 300 states are retained in every step of diagonalization. 

It is found that the low-energy state (~ 10“®iA) changes from the unscreened localized 
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Figure 6. Local density of states of the vr electrons at a point defect as a function of energy 
in the magnetic filed of (a) B/B^ua = 0, 1, and 3 and (b) = 4, 5, 7, and 10. Umin 

is determined by the system size. In the present calculation, ilmin — 21.2 T. [Reprinted from 
Taro Kanao, Doctor Thesis, Department of physics. Graduate school of science. The university 
of Tokyo, Tokyo (2014). ] 


moment state (5def//cB log 2 = 1) to the Kondo-Yosida singlet state (5def/fe_B log 2 = 0) with 
increasing hybridization from VjD = 0.01 to 0.05. Since D is about 8 eV, V is of the order 
of 0.1 eV. Namely, this Kondo effect can occur with a realistic smaller value of hybridization 
between the sp^ orbital at the defect and the vr electrons. 

In the case oiV/D = 0.1, 5def is negative in a certain range of temperature. Note that 5def is 
defined as the deference between the entropy in the presence of defect and that without defect. 
This is due to the large hybridization comparing with the band width. The similar behavior has 
been known in the narrow-band Anderson model |40j . 

In order to understand the importance of localized vr orbital in detail, we also discuss the 
width of localized vr orbital dependence of the Kondo temperature. In Fig. [5Kb), the Kondo 
temperatures Tk/D as a function of VjD for several w/D are shown. It is found that the 
Kondo temperature is sensitive to the width: With smaller w/D, Tk becomes larger owning to 
the larger LDOS at the Fermi energy. 

4. Effects of Magnetic Field 

As mentioned in the introduction, the negative magnetoresistance has been observed 
experimentally in this Kondo effect [8]. In the usual Kondo effect, the Zeeman effect on 
the localized magnetic moment causes a negative magnetoresistance, because it suppresses the 
formation of the Kondo-Yosida singlet. The Kondo-Yosida singlet is suppressed when the energy 
scales of the magnetic field and the Kondo temperature become the same order, which leads to a 
characteristic magnetic field in the negative magnetoresistance. In the Kondo effect in graphene 
with defects, however, the experimentally observed characteristic magnetic field is smaller than 
the energy scale of its Kondo temperature by one order of magnitude. This means that this 
Kondo effect (Kondo temperature) is more sensitive to the magnetic field than the usual one. 
However, the origin of this behavior has not been understood. 

In the previous section, we have shown that the Kondo screening is assisted by the peak of 














the LDOS of vr electrons which is caused by the defect scattering. Since the Kondo temperature 
sensitively depends on the height of the peak, small modulation of the localized tt orbital by 
the magnetic field will cause appreciable change in the Kondo screening. Thus, to clarify the 
modulation due to the orbital motion of the vr electrons in the magnetic field, as a first step, 
we solve the nearest-neighbor tight-binding model for the vr states of graphene in the magnetic 
field with a single point defect. In this case, the Hamiltonian is given as 

■^gra = + h.c.) , (15) 

(b> 


where ' denotes the presence of the defect, and a* and bj are annihilation operators of a vr electron 
at site i and site j on graphene, and Oij is a Peierls phase with the “string gauge” [H]. Under 
the periodic boundary condition, the possible magnetic fields are restricted to B = with 

n = 0,1,2, • • •, where Hmin is determined by the system size. In the present calculation, 61^ 
unit cells are used, which corresponds to B^i^ ~ 21.1 T. 

Figures [6] show the LDOS at the defect as a function of energy in several magnetic fields. 
We find that the height of the peak decreases and the width increases for relatively small B 
as shown in Fig. EKa). At higher magnetic helds, the peak splits into two peaks as shown in 
Fig. EKb). These results show that the localized vr orbital becomes broad and the LDOS of 
the Dirac point (the Fermi level) decreases by the orbital magnetic field. Since the Kondo 
temperature is strongly sensitive to the width of the localized vr orbital, it is expected that the 
Kondo temperature decreases drastically by the orbital magnetic field. 

5. Summary 

In conclusion, we discussed the role of localized vr orbital on the defect-induced Kondo effect in 
graphene on the basis of a numerical renormalization group study. We found that the localized 
vr orbital assists the Kondo effect due to the sp^ orbital, and the Kondo temperature is sensitive 
to the broadening of the localized vr orbital. Secondly, in order to clarify the mechanism of the 
magnetic sensitive Kondo effect, as a first step, we study the effect of the orbital magnetic field 
on the localized vr orbital by the tight-binding model with the Peierls phase. We also found that 
the LDOS of the localized vr orbital at the fermi level decreases as the magnetic field increases. 
It is suggested the the defect induced Kondo effect is strongly sensitive to the orbital magnetic 
field. 
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